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Abstract We describe how protonium, the quasi-stable antiproton-proton 
bound system, has been synthesized following the interaction of antiprotons 
with the molecular ion H^~ in a nested Penning trap environment. From a 
careful analysis of the spatial distributions of antiproton annihilation events 
in the ATHENA experiment, evidence is presented for protonium produc- 
tion with sub-eV kinetic energies in states around n = 70, with low angular 
momenta. This work provides a new 2-body system for study using laser 
spectroscopic techniques. 

Keywords Protonium • Exotic atoms • Antiprotons 
PACS 36.10-k • 34.80.Lx • 52.20.Hv 



1 Introduction 

The availability of a high-quality low energy antiproton (p) beam delivered 
by the CERN Antiproton Decelerator (AD) to the ATHENA, ATRAP and 
ASACUSA experiments has permitted the routine production of stable pure 
antimatter systems (antihydrogen, H) [TJ[2] and metastable mixed matter- 
antimatter systems (antiprotonic helium or pHe + ) [3111], both very important 
steps towards the goal of testing CPT symmetry. 

Another metastable exotic atom that is also of great interest is antipro- 
tonic hydrogen (pp), also called protonium (Pn). Its level structure is similar 
to that of hydrogen, but binding energies are much larger and its simple two- 
body nature allows an independent CPT test. In fact, spectroscopic measure- 
ments on Pn, made in near-vacuum conditions with high-precision laser tech- 
niques, would make it possible to improve the precision of the "(anti)protonic 
Rydberg" constant or, equivalently, the (anti)proton to electron mass ratio. 

Even though Pn has been studied in the past (see e.g. [5,6] ), in this paper 
we report a radically new method to produce Pn resulting in emission with 
very low kinetic energy (from some meV to ~ 1 eV) in vacuum conditions 
that could open the way to laser spectroscopic studies. Previous experiments 
produced Pn by injecting antiprotons into a molecular hydrogen target (H2), 
either liquid or gaseous. This results in a Pn lifetime which depends strongly 
upon the target density due to the effect of collisional de-excitation and which 
makes spectroscopy impossible. In the ATHENA apparatus Pn has been 
produced after a "chemical" reaction between p's and molecular hydrogen 
ions (H^ - ) trapped together with positrons (e + ) in a nested Penning trap. 
Protonium production accompanies the H production described in [T], and 
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here we describe how the two modes of p annihilation have been distinguished 
(see also [7]). 

2 Experimental details 

The ATHENA apparatus, described extensively in [8], consisted of a multi- 
electrode system of cylindrical Penning traps, 2.5 cm in diameter and ~ lm 
in length kept in an axial magnetic field of 3 T. In the 15 K cryogenic en- 
vironment of the trap, only hydrogen and helium were present in gaseous 
form, giving a residual pressure of ~ 10~ 12 Torr. Antiprotons from the AD 
were caught, cooled by electrons and stored in the so-called mixing trap. 
The latter is a nested Penning trap, approximately 10 cm long, that allowed 
e + s and p's to be confined simultaneously Under typical conditions the trap 
contained a spheroidal plasma of ~ 3.5 x 10 7 e + and ~ 10 4 p's. The resulting 
p annihilations were monitored for ~ 60 s by detectors [5] that, recording the 
passage of the charged pions, allowed reconstruction of annihilation vertices 
with an uncertainty of a few mm. 

In typical operating conditions, annhilation of p's originate from: 

— H formation followed by annihilation on the electrode surface p~ll 1 0J : 

— p annihilation in some well-defined "spots" on the electrode walls due to 
radial transport [11] : 

— annihilation following interactions with residual gas atoms or ions present 
in the trap. 

It was shown in [12] that, when the e + cloud was kept at the trap envi- 
ronment cryogenic temperature of ~ 15 K (a situation called "cold mixing" 
hereafter), annihilations were mainly due to H, even if some annihilations 
near the trap axis were present (see also FigH^). On the contrary, when the 
e + cloud was heated (by a radio-frequency drive applied to an electrode of 
the trap [I3UH]) to a temperature, T e , of several thousand K (about 8000 K 
for the data reported here, which we call "hot mixing" hereafter) H formation 
was strongly suppressed |15j and ps annihilated mainly without forming H 
(see also Fig[Tp). 

In the following, we will study in detail the latter annihilations. 

3 Results and discussion 

In FigflJ we report the x-y scatter plots (both coordinates are in the plane 
perpendicular to the trap axis) for cold and hot mixing data. In Fig[2l the cor- 
responding r-z scatter plots are reported (where r — ^/ x 1 + y 2 is the radial 
position, i.e. the distance from the trap axis, and z is the axial coordinate, 
measured from the symmetry plane of the trap). 

Even though all the distributions are broadened by the uncertainty in the 
vertex reconstruction, it is clear that two different structures are merged in 
cold mixing case, while in hot mixing just one appears. For the cold mixing, 
besides the annihilations on the trap wall situated at r = 1.25 cm due mainly 
to H (as shown in |12j ) and having a relatively wider z-distribution (10j . 
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Fig. 1 Scatter plot of annihilation vertices projected on a plane perpendicular to 
the trap axis for (a) cold mixing, (b) hot mixing. The dashed line represents the 
trap wall. 



there are some annihilations situated at smaller r and with a very sharp 
z-distribution (see also FigO)). For hot mixing only the latter are present, 
though their axial distribution is broader, as is the radial distribution (see 
also FigHk,b). 

The capability of the ATHENA detector to detect the spatial and tem- 
poral coincidence between p and e + annihilations, and therefore to separate 
H from other annihilations, allows us to infer that for hot mixing almost all 
annihilations (even near the wall) are not due to H, while for cold mixing 
we have to distinguish between annihilations near the trap axis, that are not 
due to fi (apart from a few poorly reconstructed vertices), and near the wall, 
where the number of non-fi annihilations is negligible. 

This is clear if we look at FigJ31 where we consider annihilations happen- 
ing in coincidence with two (and only two) photons detected, and we plot 
the cosine of the angle between the two detected photons, cos(# 77 ). This 
distribution should have a peak in cos(# 77 ) = — 1 if annihilations are due to 
H, because of the two back-to-back 511 keV photons produced by the e + 
annihilation, while we do not expect any peaks for annihilations not related 
to H. The former is the case for cold mixing on the wall (Fig(3^,), whilst 
the latter is the case for cold mixing near the trap axis (Figj3]b) and for hot 
mixing (Fig|3b). 

It is notable that the axial extent of the annihilation distributions is 
small compared with the size of the nested trap, so they cannot be in-flight 
annihilations of p's on residual gas. However, their radial extent is too large 
to be explained by in-flight annihilations on positive ions trapped inside the 
positron well. So, only a process involving p interaction with an ion to form 
a neutral system that can annihilate in flight, perhaps after leaving the trap, 
can explain the data. 

Time-of-flight measurements following charged particle ejection from the 
trap excluded the presence of protons and of ions of atoms more massive than 
helium. However, reactions of p with helium ions would give distributions that 
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Fig. 2 Scatter plot of annihilation vertices as a function of r and z for (a) cold 
mixing, (b) hot mixing. The dashed line represents the trap wall. The semi-ellipse 
indicates the size and position of the positron cloud. 
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Fig. 3 Histogram of the cosine of the angle between the two detected photons 
(cos(# 77 ); see text) for (a) cold mixing, on the wall (r > 1 cm); (b) cold mixing, 
near the trap axis (r < 0.5 cm); (c) hot mixing. 



do not match the data, since p He ++ gives rise to a charged system followed 
by rapid annihilation, and also for p He + the residual electron would be 
ejected in less than 10 ns |16) . again giving rise to a charged system. 

Further important information on the p-system formed can be found by 
exploiting the number of tracks from each annhiliation vertex, corresponding 
to the number of charged pions produced during the annihilation, since this 
depends on the particle with which the p annihilates (see e.g. [17] ). TabQ] 
shows the ratios, R23, of the number of the reconstructed annihilation ver- 
tices having two tracks to those with three tracks, for different data samples. 
In order to have a better understanding of these data, a Monte Carlo simu- 
lation of pp (i.e. Pn) annihilations inside the ATHENA apparatus has been 
performed. Comparing these data, we see that annihilations on wall differ 
from those originating inside the trap, which are compatible with pp annihi- 
lations both for cold mixing and hot mixing. 



6 



Table 1 Experimental and Monte Carlo results for the number of charged pion 
tracks due to p annihilations. 



Data set 


Ratio R23 


on wall Ratio R23 at centre 


Cold mixing 


1.35±0.01 


1.22±0.04 


Hot mixing 


1.38±0.10 


1.17±0.04 


ps only (no mixing) 


1.40±0.03 




Monte Carlo pp 




1.19±0.01 



Combining this information, we infer that the most probable p-ion reac- 
tion is: 

p + Hj -> Pn(n,0 + H. (1) 

where the Pn leaves the trap (because it is neutral) , and has an exponentially 
distributed annihilation lifetime (depending on n and I). For relatively high 
n, Pn is formed in a metastable state and its lifetime can be greater than 
l//s Q2]. 

The ions required for reaction ([T]) may have been created during the 
positron loading procedure [19] . due to collisions with the residual H2 gas, or 
during p loading (similarly to that reported in [20] ). Measurements of charge 
during the dump of the trapped particles has indicated that the number of 
ions could be as high as 10 5 , depending on the vacuum conditions. 

In order to check our hypothesis, we performed a self-standing Monte 
Carlo simulation to reproduce the observed annihilation distributions for hot 
mixing and cold mixing in the following way: 

— we used the information we have on the e + plasma shape |13lH"4] to gen- 
erate the Pn starting point distributions (in particular, it was found that 
the e + plasma was approximately a spheroid with radius r p =l mm and 
axial half-length z p =\6 mm, rotating with a frequency of 300 kHz; i.e. a 
surface velocity of about 2000 ins -1 ). 

— assuming that Pn is produced in thermal equilibrium with the e + plasma, 
and knowing that Pn must inherit the drift velocity from its charged 
components, we generated a velocity distribution summing the thermal 
Maxwellian isotropic velocity (with a mean value fixed by the measured 
temperature of the e + plasma) with a velocity along the tangential direc- 
tion (that is the same as inferred by the measured parameters of the e + 
plasma, since the drift doesn't depend on either the mass or the charge 
of the particle). 

Starting with hot mixing (8000 K), the simulation was performed assum- 
ing that the Pn was produced on the surface of the e + spheroid (character- 
ized by the parameters above), with a Gaussian distribution along z (having 
a = 10 mm and being limited to \z\ < z p ) or alternatively inside the e + 
spheroid with a uniform density. The result is nearly independent of the as- 
sumed distribution of starting positions, since the system is dominated by the 
thermal velocity of 5600 ms . The simulated radial and axial distributions 
for the best fitted mean lifetime (1.1/xs) are superimposed on the experimen- 
tal data in FigHJ The agreement is good, and in particular the simulation 
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Fig. 4 Distributions of annihilation vertices for hot mixing (continuous line) , with 
superimposed Monte Carlo simulation results (dashed line), (a) radial distribu- 
tion; (b) axial distribution, for events near the trap axis (r < 0.5 cm); (c) axial 
distribution, for events near the trap wall (r > 1 cm). 



shows that about 25% of the Pn produced reaches the wall, as observed, and 
it also predicts that the axial distribution near the wall is slightly wider than 
that near the axis (Fig@fc>,c) due to the essentially isotropic spreading. 

For cold mixing the analysis is less straightforward, because the H con- 
tribution on the trap wall must be subtracted. To do so, we considered the 
difference between a radial distribution taken in a z-slice where Pn is present 
(e.g. \z\ < 0.5), and one where there is no Pn (e.g. 0.5 cm < \z\ < 1.5 cm), 
normalized on the tail for r > 1.5 cm, which is essentially H only. The re- 
sult is shown in Fig[5^,. For the axial distribution, isolating the Pn signal is 
much simpler, because only annihilations near the axis of the trap need be 
considered (Fig(5b). 

In order to simulate the cold mixing distributions, the parameters (r p , z p , 
frequency rotation) of the e + plasma were the same as for hot mixing, and 
we generated Pn with a thermal velocity corresponding to 15 K (250 ms" 1 
along each direction). However, in this case it is necessary to assume that Pn 
is produced in a very narrow region around r = r p and z = (we generated 
it with a Gaussian distribution along z with a = 2.5 mm). The best fitted 
mean lifetime is 1.1 /is, as found for hot mixing. 

The agreement between the experimental and Monte Carlo data is good 
(Figj5^,b). Furthermore, Monte Carlo results show that less than 0.5% of 
Pn reached the wall in this case, thus confirming the consistency of our 
normalization technique. 

For cold mixing, we were also able to analyze another data sample, where 
the positron plasma had different parameters (r p =2.5 mm, z p =18 mm, rota- 
tion frequency of 80 kHz, surface velocity of about 1300 ms _1 ). Following the 
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Fig. 5 Distributions of annihilation vertices for cold mixing, (continuous line), with 
superimposed Monte Carlo simulation results (dashed line), (a) radial distribution, 
after removing the H contribution as discussed in the text; (b) axial distribution, 
for events near the trap axis (r < 0.5 cm); (c) radial distribution, for the cold 
mixing sample with wider positron plasma (in this case, the thick dashed line 
corresponds to Monte Carlo results with the correct plasma parameters, while the 
thin dashed line corresponds to Monte Carlo results with the parameters of the 
narrower positron plasma; see text). 



same prescription, we have reproduced the radial distribution of this sample 
(Fig{5fc), showing there is a notable difference between the distribution of Pn 
annihilations in this case and that for a narrower e + plasma. 

It has to be stressed that it is impossible to reproduce the experimental 
results unless Pn is produced in a narrow region straddling the "equator" of 
the e + spheroid. A possible and quite straightforward way to explain this, is 
to assume that there is some kind of (partial) separation between the e + and 
the Hj" ions, as expected in case of thermal equilibrium (sec |21,22J), and as 
has been observed for a mixture of e + and 9 Be + in 23 , 24J . In fact, for our 
experimental conditions, assuming thermal equilibrium, HJ" ions experience 
a centrifugal potential barrier of the order of 10 meV meaning that, for 15 K, 
their thermal energy (~ 1 meV) is not enough to allow the ions to penetrate 
the e + plasma. However, for 8000 K, the barrier is negligible compared to the 
thermal energy (~ 700 meV) of the ions such that they will be distributed 
uniformly inside the plasma. 

Another constraint from our Monte Carlo simulation relates to the Pn 
kinetic energy. In fact, the hot mixing and the cold mixing data cannot be 
reproduced using the same mean lifetime for both of them if the system has 
a recoil kinetic energy of the order of 1 eV or greater. This probably means 
that the p-Hj collision happens preferentially as a resonant transfer, i.e. the 
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dissociation energy of Hj" is taken care of by the binding energy of Pn, so 
that its energy level, n, should be around 70. 

Because we find a mean lifetime of about 1.1/is, it can be inferred [18] 
that the orbital angular momentum, Z, should be around 10. This could be a 
consequence of the fact that, with such a slow relative collisional velocity, the 

molecular ion will be strongly polarized, giving rise to an almost collinear 
collision. 



4 Conclusions 

In this paper, we have presented evidence for the production of protonium 
in vacuum. It is formed in a metastable state (with a lifetime of about 1.1/xs, 
independent of the environment temperature) and with near-thermal kinetic 
energies (varying from some meV to less than 1 eV). The number of pro- 
duced protonium atoms was about 100 for each mixing cycle (~ 100 s), in 
which around 10 4 p's were injected into the mixing trap, while the estimated 
number of ions trapped with the positrons was typically 10 4 - 10 5 . Taking 
into account the recently achieved capability of accumulating ~ 10 8 and 
storing ~ 5 xlO 6 p's in some minutes [25inZo] . our result opens up the possi- 
bility of performing detailed spectroscopic measurements on protonium as a 
probe of fundamental constants and symmetries. 
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